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Abstract. Two neutrino double beta decay of ' 6 Ge to excited states of ,6 Se has been 
studied using data from Phase I of the Gerda experiment. An array composed of up 
to 14 germanium detectors including detectors that have been isotopically enriched 
in 76 Ge was deployed in liquid argon. The analysis of various possible transitions to 
excited final states is based on coincidence events between pairs of detectors where a 
de-excitation 7 ray is detected in one detector and the two electrons in the other. 

No signal has been observed and an event counting profile likelihood analysis has 
been used to determine Frequentist 90 % C.L. bounds for three transitions: 0+ s — 2f: 
Tfj 2 > 1.6-10 23 yr, 0+ s - 0+: T?f 2 > 3.7-10 23 yr and 0+ s . - 2+: Tfj 2 > 2.3-10 23 yr. 
These bounds are more than two orders of magnitude larger than those reported 
previously. Bayesian 90 % credibility bounds were extracted and used to exclude several 
models for the 0 + s — 0 ^ transition. 


PACS numbers: 23.40.-s, 21.10.Tg, 27.50.+e, 29.40.Wk 
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1. Introduction 

The observation of neutrinoless double beta (0 vfdfd) decay would imply physics 
beyond the Standard Model of particle physics because the process manifests lepton 
number violation. Assuming light neutrino exchange as the dominant process, the 
experimentally observable half-life would be connected to the effective Majorana 
neutrino mass via a phase space factor (F) and a nuclear matrix element (NME, M). 
Calculations of M. depend strongly on the nuclear structure model chosen. Uncertainties 
in the model translate to uncertainties when converting a measured half-life into the 
Majorana neutrino mass or its limit respectively [1]. 

Two neutrino double beta (2 u/3/3) decay is a weak second order Standard Model 
process which has been observed in eleven nuclides and double neutrino electron capture 
(2z/ECEC) in two nuclides with half-lives between T 1 2 ^ 2 =10 18 — 10 24 yr [2, 3]. Due to the 
different reaction mechanisms the effective Majorana neutrino mass (mpp) enters only 
in the half-life for the (V/3/3 mode, and not for the 2 v(3(3 mode: 

2u/3/3 : (T 2 /^)' 1 = F 2v ■ \M 2u \ 2 (1) 

Oz/Ad : C Z?; 2 ) -1 = F 0u • \M 0v \ 2 • (myp) 2 (2) 

The NME for 2 u/3/3 and 0 u/3/3 (M 2u , M Qv ) and the respective phase space 
factors (F 2u , F 0u ) are numerically different but rely on similar model assumptions. 
Experimentally verifying the calculations of Eq. (1) reduces to some extend uncertainties 
of calculations of A4°U For a recent review see Ref. [4], 

Apart from decays into the ground state (g.s.), double beta decays can also occur 
into excited states of the daughter nucleus. These decay modes are expected to have 
a lower rate due to a smaller phase space. However, their experimental signature is 
enhanced by the accompanying de-excitation 7 rays. Excited state transitions can 
in principle occur in both the 2 u/3/3 and the 0 u/3/3 mode distinguishable only in the 
different sum of residual electron energies. These investigations provide additional 
information on the nuclear structure models and more experimental constraints on the 
system of Eqs. (1 and 2 ). So far only transitions to the first excited 0+ state have been 
observed, first for 100 Mo [5] in 1995 and for 150 Nd [ 6 ] in 2004. Recent half-life values are 
T*j 2 = (7.5 ± 1 . 2 )- 10 20 yr [7] and T*f 2 = (1.33+^)-10 20 yr [ 8 ], respectively. 

The double beta decay of ,6 Ge to the ground state of 76 Se has a Q-value of 
Q/3/3=(2039.061±0.007) keV [9] and can potentially feed any 0 + or 2 + excited state in 
76 Se up to this energy. The search in this work focuses on the 2 u/3/3 transitions from the 
ground state 0+ s of 76 Ge to the three lowest excited states in 76 Se: 2f, Oj 1 ", and 2 2 (see 
Fig. 1). The phase spaces of 2 u/3/3 transitions scale with the total available energy E 
as F 2u ~ E 11 and reduce the rates for higher energetic states. The rate is further 
suppressed by the spin-constraint for 2 + states. For the investigated transitions the 
existing experimental upper limits for half-lives are shown in Table 1. In addition results 
from theoretical calculations based on various nuclear structure models are presented. 
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Figure 1. Possible double beta decay modes of 76 Ge to excited states of 76 Se that 
are investigated in this work. 


The largest rate is expected for the 0^ s — 0^ transition and lies within the experimental 
sensitivity of this analysis. Three calculations have been recently performed for this 
transition. One is based on the renormalized proton-neutron QRPA using wave functions 
from Ref. [10]. This calculation predicts a half-life of Tjy 2 (0+)=(1.2 — 5.8) • 10 23 yr for 
axial vector couplings of g A =1.2Q - 1.00 [11]. Employing a microscopic interacting 
boson model (IBM-2), the half-life ratio between the 0+ s — 0^ and the ground state 
transition 0+ s — 0+ s is calculated with NMEs and phase space factors from Refs. [12, 13]. 
The predicted half-life ratio is 3300 for 76 Ge [14]. Scaling this ratio with our recently 
measured ground state half-life of T 1 2 / l/ 2 (0+ s ) = (1.926 ± 0.095) • 10 21 yr [15] results in the 
predicted excited state half-life T^ 2 (0f)=6A ■ 10 24 yr independent of g A . A shell model 
(SM) calculation predicts 3—6.7)-10 24 yr [16] assuming g A =1.2Q - 1.00. The 

range of Tff 2 (Qf) additionally encompasses results of two different effective interactions 
that have been used also in Refs. [17] and [18]. The SM calculation of TfrM) to 
the ground state is in good agreement with the experimental data. The IBM-2 and 
SM predictions are significantly longer than other calculations. The current status for 
the experimental and theoretical situation is summarized in Table 1 for all investigated 
transitions. 

2. The coincidence analysis of Gerda Phase I data 

The 7 cascade following excited state transitions provides a well defined experimental 
signature which enables large background suppression. The granular installation of the 
Gerda setup is used to measure coincidences between two germanium detectors. In 
the following the two detectors are distinguished between ( 1 ) a “source” detector where 
the 2 u/3/3 decay occurs and the two electrons are detected and ( 2 ) a “gamma” detector 
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Table 1. Experimental half-life limits are compared to model predictions for 2v/3/3 
excited state decay modes in 76 Ge as discussed in this paper. The energy of the final 
level is given. Abbreviations denote: HFB: Hartree-Fock-Bogoliubov, QRPA: Quasi 
Random Phase Approximation, MCM-QRPA: Multiple Commutator Model QRPA, 
RQRPA: Renormalized proton-neutron QRPA, IBM: Interactive Boson Model, SM: 
Shell Model. 


2i//3/3 decay mode 

Tfj 2 [yr] 

mo del/exp. 

ref. 

year 

0+ s . - 2+ (559.1 keV) 

> 6.3-10 20 (68 % C.L.) 

exp. 

[19] 

1992 


> 1.1-10 21 (90 % C.L.) 

exp. 

[20] 

1995 


1 .2-10 30 

SM 

[21] 

1984 


5.8-10 23 

HFB 

[22] 

1994 


5.0-10 26 

QRPA 

[23] 

1994 


2.4-10 24 

QRPA 

[24] 

1996 


7.8-10 25 

MCM-QRPA 

[25] 

1996 


1 .0-10 26 

RQRPA 

[26] 

1997 


(2.4 — 4.3)-10 26 

RQRPA 

[27] 

1998 


2 .0-10 27 

RQRPA 

[28] 

2014 

0+ s - 0+ (1122.3 keV) 

> 6.3-10 20 (68 % C.L.) 

exp. 

[19] 

1992 

> 1.7-10 21 (90 % C.L.) 

exp. 

[20] 

1995 


> 6.2-10 21 (90 % C.L.) 

exp. 

[29] 

2000 


1.32-10 21 

HFB 

[22] 

1994 


4.0-10 22 

QRPA 

[23] 

1994 


4.5-10 22 

QRPA 

[24] 

1996 


7.5-10 21 

MCM-QRPA 

[25] 

1996 


(1.0- 3.1)-10 23 

RQRPA 

[26] 

1997 


(1.2 — 5.8)-10 23 

RQRPA 

[11] 

2014 


6.4-10 24 

IBM-2 

[14, 15] 

2014 


(2.3- 6.7)-10 24 

SM 

[16] 

2014 

0+ s , - 2+ (1216.1 keV) 

> 1.4-10 21 (90 % C.L.) 

exp. 

[20] 

1995 


1 .0-10 29 

QRPA 

[23] 

1994 


1.3-10 29 

MCM-QRPA 

[25] 

1996 


(0.7- 2.2)-10 28 

RQRPA 

[26] 

1997 


where the de-excitation 7 ray is detected. This is achieved by searching for a 7 ray of 
interest in one detector and labeling the other one as “source”. Note, that the distiction 
is made on the analysis level and that it might not be unique in some rare cases. 

2.1. Event signature 

The signatures of the investigated transitions are listed in the following (see also Fig. 1). 

(i) The transition feeding the 559.1 keV level (0+ s — 2+) has one single de-excitation 
7 ray of the same energy and thus a f3f3 spectrum with 1480.0 keV endpoint energy. 

(ii) The transition feeding the 1122.3 keV level (0+ s — 0] 1 ") de-excites via the 2\ 
state. A 563.2 keV 7 ray is followed practically immediately by the 559.1 keV 7 ray. 
The angular correlation W of the two 7 rays in the 0^ — 2\ — 0+ s cascade is given by 
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W(9) oc 1 — 3 cos 2 9 + 4 cos 4 6 where 9 is the angle between them [30] . The analysis 
searches for one of the two 7 lines in the “gamma” detector. However, the energy 
resolution in Gerda does not allow to separate the two 7 lines and a single peak region 
is used as region of interest. The (3/3 endpoint energy is reduced to 916.8 keV. 

(Hi) The transition to 1216.1 keV level (0^ s — 2f) has two decay branches: branch 1 
with two de-excitation 7 rays via the 2f state with energies of 657.0 keV and 559.1 keV. 
Branch 2 has a single de-excitation 7 ray of 1216.1 keV directly into the ground 
state. The former has a branching ratio of 64 % leading to 36 % for branch 2. The 
angular correlation W of the two 7 rays in the 2\ — 2f — 0+ s cascade of branch 1 is 
W(9) oc 1 — j|cos 2 0 + ]|cos 4 0. The (3(3 endpoint is reduced to 822.0 keV. The two 
branches are treated separately in the analysis and are later combined into a single value 
for the half-life. 

2.2. The Gerda experiment 

The GERmanium Detector Array (Gerda) is an experiment designed to investigate 
neutrinoless double beta decay of 76 Ge with an array of high purity germanium (HPGe) 
detectors made from material with the ,6 Ge fraction enriched to ~87%. The detectors 
are operated within a cryostat containing 64 m 3 of liquid argon (LAr) at the Laboratori 
Nazionali del Gran Sasso (LNGS) of Istituto Nazionale di Fisica Nucleare (INFN). The 
installation of the HPGe detectors in a closely spaced array with little material between 
the detectors facilitates an anti-coincidence veto for the 0 v(3(3 search into the ground 
state suppressing 7 ray background. 

The Gerda setup is described in detail in Ref. [31]. The array of HPGe detectors 
used in Phase I of the experiment was organized into 4 strings of 3 to 5 detectors each. 
Three different detector types were used: semi-coaxial and BEGe detectors enriched 
in 76 Ge ( enr Ge) and semi-coaxial detectors with natural abundance ( nat Ge). Five enr Ge 
semi-coaxial detectors are from the Heidelberg-Moscow experiment [32], three enr Ge 
semi-coaxial detectors from the Igex experiment [33], and three nat Ge semi-coaxial 
detectors from the GENIUS test facility [34], All were refurbished for their operation in 
LAr. Additionally, an initial batch of five enr Ge BEGe detectors produced for Gerda 
were employed [35]. The detector strings were lowered into LAr by a two-arm lock 
system supporting one string in one arm and three strings in the other [31]. 

2.3. Data sets and energy calibrations 

The data of Gerda Phase I are used in this analysis. As described in Ref. [36], some 
detectors were removed and replaced by the 5 BEGe detectors after an initial data 
taking period. A higher background was observed for 49 d due to this change (silver 
data set as defined in Ref. [37]). For the present analysis this higher background period 
is omitted. The periods before and after had different detector array configurations and 
are treated independently in the analysis. 
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The data set is composed of runs of approximately 1 month each. Each selected run 
shows a stable operation of the whole array [37]. In some runs, individual detectors are 
excluded from the analysis due to temporary instabilities. In the present analysis every 
detector - also if its data is not used - is considered as source of the decay. A two-detector 
coincidence may be registered also in case the decay occurs in an excluded detector and 
the two 7 rays are deposited in two active detectors. The efficiency calculation is taking 
these possibilities properly into account. 

In the analysis, the target mass is defined as constant within each array 
configuration; changes due to inactive detectors are condensed into the signal detection 
efficiency of the array. The live-time sums up to an exposure of £=31.04 kg-yr for the 
whole data set including all nat Ge and enr Ge detectors. The isotopic exposure of ,6 Ge is 
£76=22.3 kg-yr. 

The same energy software-threshold of 100 keV is applied to all detectors. This 
basic threshold was applied to ensure a full reconstruction efficiency for all detector 
energies. After additional quality cuts and a /i-veto cut the data set contains 2710 
two-detector events and 82 three-detector events. This can be compared with «7 ■ 10 5 
single-detector events. The efficiency to detect de-excitation 7 rays in three-detector 
events is more than one order of magnitude smaller and thus only two-detector events 
are further analyzed. 

Energy calibrations were performed with 228 Th sources typically once per week. 
In case of an energy deposition in more than one detector, cross-talk affects the 
reconstructed energy. Data taken by dedicated calibrations and the 42 K 7 line in 
the physics data allow to measure the effect and hence to correct for it. The energy 
dependence of the cross-talk is linear in good approximation. The exposure averaged 
energy resolution (FWHM) of the 583 keV 7 line from 208 T1 is 4.2 keV for coincident 
events while it is 3.8 keV for events with energy deposition in one detector only. The 
uncertainty on the resolution is estimated to be 10 %. 

2-4- Monte Carlo simulations 

Monte Carlo (MC) simulations were used to construct a specific background model 
for coincidence events and to determine the detection efficiencies. The background 
model described in Ref. [37] is used as a starting point whereto additional background 
contributions had been added. Due to the coincidence requirement the individual 
background contributions to the spectra are of different significance when compared 
to the single-detector spectrum. Therefore, the evaluation of the background sources 
was tuned specifically for two-detector events. The main sources of background are 
contributions from 214 Bi, 212 Bi, 208 T1, 228 Ac, 40 K and 60 Co on the detector holders, 42 K 
and 39 Ar distributed homogeneously in the LAr, 42 K on the detector n + contact and 214 Bi 
on the detector p + contact. Additionally, one further background contribution from 
108mAg j n t} ie signal cables was considered which is not part of the minimal background 
model in Ref. [37] and only visible in coincidence data. 
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The energy spectra from the data obtained with all the detectors and from the 
background model are shown in Fig. 2 for two-detector events. This background 
description is used for the development of cuts and the calculation of sensitivities but it 
does not enter itself in the final analysis. The left panel of Fig. 2 shows the individual 
detector energies in which each of the two detectors has a separate entry. The right panel 
shows the sum-energy spectrum; i.e., the total energy deposited in the array for a two- 
detector event. The shaded histograms show the data which can be directly compared 
to the background model shown in black. The individual background components for 
the semi-coaxial detectors only are shown in color; due to low statistics the ones for the 
BEGe detectors are omitted in the plot for clarity. 
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Figure 2. Data and background model for two-detector events after an individual 
detector threshold of 100 keV before cut optimization. In the single-energy spectra 
(left) each detector has a separate entry per event in the histogram whereas in the 
sum-energy spectra (right) the two energies are summed before histogramming. Data 
events are shown in solid gray and the background components for the coaxial detectors 
by the colored lines. 


There are some differences found between background MG and data. 

At low energy beta decay events from 39 Ar with end point energy 565 keV dominate. 
The detected energy is highly sensitive to the exact knowledge of the detector dead layer 
which is not available. Hence data and Monte Carlo disagree here. The probability of a 
two-detector coincidence due to 39 Ar is small and highly sensitive to the exact thickness 
of the detectors’ dead layer. However, this effect is not relevant for the energy regions 
investigated in this analysis. Another excess is visible around 1.8 MeV in the single¬ 
energy spectrum in the data. However, it can not be due to a missing 7 line or beta 
spectrum in the background model and furthermore the exess is not very significant. 
The agreement between the model and the data is sufficient for this analysis since the 
former is only used for cut optimization [38]. 
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MC simulations were used to determine the signal efficiency for each decay mode. 
The simulations were performed with MaGe [39] which considers the angular correlation 
between the de-excitation 7 rays. The detectors that were inactive in a given run are 
also set as inactive in the post-processing of the MC data. Each detector pair in each 
run has an individual signal detection efficiency which is then life-time weighted into a 
single number for the whole data set. 

3. Analysis 

The analysis is based on event counting in a region of interest (ROI) after a sequence 
of cuts. This is not a blind analysis, but special care was taken to avoid biases. The 
construction of the sequence of cuts starts with initial choices that are then optimized 
to maximize sensitivity. This aims to prevent ad hoc choices of analysis parameters. 
The MC background model and efficiency are used for the cut optimization. Systematic 
uncertainties and potential deficiencies in the background model only affect the choice 
of analysis parameters and do not have a direct effect on the derived half-life results 
since the background in the ROI is estimated from side bands (SB). 

The spectra of the simulated 2 u/3/3 decays scaled to 10 23 yr half-life are shown 
in Fig. 3. Also shown is the corresponding background model. The single-energy 
spectra and the sum-energy spectra are shown separately. The 2 u/3/3 decays have a 
continuous shape in the sum-energy spectra due to the continuous electron component 
that is almost always detected in the source detector. The single-energy spectra show 
distinct de-excitation 7 lines being detected outside the source detector enabling strong 
background reduction. 



energy [keV] 


energy [keV] 


Figure 3. Illustration of the 2v/3/3 decays for two-detector events scaled to a half-life 
of 10 23 yr for each decay mode. Shown are the single-energy spectra (left) and the 
sum-energy spectra (right) befor cut optimization. Also shown are the background 
model (black line) and data (gray histogram). 
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The energy distributions of two-detector coincidence events are shown in scatter 
plots (Fig. 4) for the simulated 0+ s — Oj 1 " transition, the background model and the 
data. The color scale denotes the expected event densities for an excited state half- 
life of 10 23 yr and for the expected number of background events in the data set, 
respectively. The black points are data events. In this representation, horizontal and 
vertical lines are features of the single-energy spectra whereas diagonal lines are features 
of the sum-energy spectra. Many high energy background 7 lines originate from outside 
the germanium detectors and can scatter into two detectors. These events can be 
suppressed with a cut on the sum energy. In Fig. 4 the diagonal 7 line from 42 K is 
clearly visible in the data and the background model. Additional 7 lines from 208 T1, 40 K 
and 211 Bi are only visible in the projections (see Figs. 2). 
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Figure 4. Scatter plot for 
coincident events showing the 
simulated 2n/3/3 signal process 
for the 0+ s — Oj 8 decay mode 
scaled to 10 23 yr half-life (top 
left), the simulated background 
model (top right) and the data 
events (bottom left). The 
number of MC events is scaled 
to the Phase I data set. 


3.1. Sequence of cuts 
A sequence of four cuts is applied: 

1. standard cuts: quality cuts and // veto as in Gerda Phase I 

2 . coincidence cuts: specific two-detector cuts optimized for each decay mode 

3. background cuts: exclusion of background 7 lines 
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4. detector pair cuts: cuts to select detector pairs that exhibit a high efficiency 

The cuts are designed to optimize the sensitivity Ti/ 2 , expressed by the figure of 
merit S = ej\/~B with the signal detection efficiency e and the number of surviving 
background events B. The efficiency and the background in the sensitivity study are 
entirely based on MC simulations [40]. 

The standard cuts include quality cuts and the muon veto cut and are initially 
applied to the data set. They are identical to the cuts in the Gerda 0 u/3/3 analysis 
(without pulse shape cuts) which are described in Ref. [36] and references therein. 

The coincidence cuts are specific for each decay mode. They require a coincident 
event with either of the two detectors having the full energy of any de-excitation 7 
ray within a peak (energy) window size (PWS) of ±2 <r E where oe denotes the energy 
resolution at energy E. The PWS defines the specific ROI and SBs. Furthermore, a 
sum energy limit of 2039 keV is applied to exclude events with a total energy deposition 
larger than the possible energy release in , 6 Ge. 

SBs are created by shifting the energy window of the coincidence cuts to lower and 
higher energies. The window size of a single SB is the same as for the corresponding 
ROI. To reduce the statistical uncertainty in the background estimation, a total of 4 
SBs are defined for each decay mode; two at lower energies and two at higher energies. 
Based on the MC background model the SB regions are chosen for each decay mode 
individually. They are as close as possible to the ROI and avoid background 7 lines in 
the single-energy spectra. 

For two-detector events it is possible that the same event is tagged in the ROI and 
in a SB. A signal event that is tagged in the ROI has a 4% chance to be tagged also 
in one of the SBs. A background event in a SB has a 2 % chance to be tagged also in 
another SB or the ROI. Events that are tagged in more than one of these five regions 
are rejected to avoid double counting in the statistical analysis. This reduces the signal 
efficiency by 4 % and the background by 2 % [40]. One event in the ROI of 0^ s — is 
removed by this requirement. 

The background cuts are motivated by background 7 lines that interfere with 
the ROIs. Fig. 5 shows the simulated sum-energy spectrum of the background model 
and the signal process after the coincidence cut for the 0+ s — 0^ transition. The low 
side cutoff is created by the smallest 7 ray energy of interest (559.1 keV) in one detector 
and the single detector threshold of 100 keV in the other. The high cutoff is given by 
the sum energy limit. The dips in the spectra are created by rejecting events that are 
tagged in the ROI and a SB or two SBs simultaneously. The peak at 1122.3 keV for the 
signal process (green) is created by decays that occur in the dead layer of a detector or 
inside an excluded detector. In this case the /3 energy of the event is not detected and 
the two 7 rays trigger a two-detector event with discrete sum energy. 

The strongest peaks in the background spectrum (black line in Fig. 5) are found at 
1524.7 keV and 1460.8 keV belonging to 42 K and 40 K, respectively. They represent the 
region in the scatter plot (Fig. 4, bottom) where the diagonal sum energy 7 lines cross the 
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horizontal and vertical single energy lines of the 2 v/3/3 decay. The next strongest feature 
in the background model is the peak around 1170 keV which is a combined structure: 
Firstly, a 609.3 keV 21i Bi 7 ray coincides with another 211 Bi 7 ray that scatters and 
deposits energy around the ROI of 559.1 keV in one detector; the 609.3 keV 7 ray is 
fully detected in the other detector creating a sum energy around 1170 keV. Secondly, 
a similar coincidence with a 614.3 keV 108m Ag 7 ray and another one from the 108m Ag 
decay can happen. The other peaks originate mainly from 211 Bi, 228 Ac or 60 Co. They 
are created either ( 1 ) by a single 7 ray scattering into two detectors or ( 2 ) by two 7 rays 
with one 7 ray fully detected in one detector. In both cases an energy around ~ 560 keV 
has to be deposited in at least one detector. This can be identified in Fig. 4 where the 
horizontal/vertical cut window crosses anti-diagonal background lines in scenario ( 1 ) or 
horizontal/vertical background lines in scenario ( 2 ). 



Figure 5. Simulated sum-energy spectrum after coincidence cuts for the 0+ s — 0^ 
transition. The background model (black) and signal process scaled to 10 23 yr half-life 
(green) are shown. 

The background cuts are applied to all decay modes in the same way. 42 K 
and 40 K are excluded by requiring the detector sum energy not to be in the range 
of 1524.7 ± 5.5 keV nor in 1460.8 ± 5.4 keV. The energy ranges correspond to 2<7e- 
Additionally, the individual detector energy is required to be outside 611.0 ±5.7 keV to 
exclude the 214 Bi and 108m Ag background. 

At this stage, the combination of coincidence cuts and background cuts can be 
optimized in three ways. The individual detector threshold is scanned from 100 keV 
to 500 keV in steps of 50 keV. This is equivalent of increasing the low energy cutoff in 
Fig. 5. The figure of merit S is determined for each step and the individual detector 
threshold with the highest S is chosen as the optimized cut. In a similar way the sum 
energy limit is scanned from 1000 keV to Qpp in steps of 50 keV. This is equivalent of 
changing the high energy cutoff in Fig. 5. Finally also the peak window size is scanned 
from la to 2a energy resolution in steps of 0.2 keV. The optimization of these cut 
parameters is performed individually for all three decay modes. The optimal values are 
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shown in Table 2 along with their relative sensitivity gain compared to the base values. 


Table 2. Optimized values for individual detector threshold (IDT), sum energy limit 
(SEL) and the peak window size (PWS). The respective ROIs are given as absolute 
ranges; the * indicates an expanded range due the occurrence of more than one 7 line. 
The sensitivity gain is related to the base values before background cuts. 


mode 

IDT 

[keV] 

SEL 

[keV] 

PWS 

[keV] 

size ROI 
[keV] 

sensitivity gain 
[%] 

07 - 2 1 

450 

1750 

± 2.6 

[556.5 - 561.7] 

41.3 

o+ s . - 0 + 

250 

1750 

± 2.2 

[556.9 - 565.4]* 

18.9 

0 +,. - 2 + B 1 

250 

1800 

±2.4 

[556.7- 561.5] 

20.4 




±3.0 

V[654.0 - 660.0]* 


OJ,, - 2+ B2 

300 

1850 

± 2.8 

[1213.3 - 1218.9] 

98.5 

base values 

100 

2039 

2 <7 E 




The detector pairs of coincidence events which are used in the analysis are also 
tuned to maximize S. Each detector pair is composed of the gamma detector and the 
source detector which are identified by containing the energy deposition in the ROI 
or an arbitrary energy deposition within the allowed energy range, respectively. This 
results in an asymmetric effect for nat Ge and enr Ge detectors: enr Ge detectors are more 
likely to be the source of a 2 u/3/3 decay than nat Ge detectors but have an equal chance 
to detect 7 rays. Taking into account all data sets there are in total 28 detector pairs in 
the Gerda array and 56 detector pairs if the distinction between source detector and 
gamma detector is taken into account. The maximization of S selects 37 pairs. The 
sensitivity gain from the pair selection is between 7 and 10 % depending on the decay 
mode. More detailed information on the sequence of cuts are reported in Ref. [40]. 

After applying the full sequence of cuts, the validity of the SB regions is tested. For 
these regions a flat background is required. In case the SBs are symmetrically placed 
around the ROI, also a linear background is sufficient. Table 3 shows the expected 
and the observed event count for each SB. For the 2f mode the model predicts a flat 
background while for the others the expected counts at higher energies (SB4) is slightly 
smaller. However the average agrees with the background model prediction in the ROI. 
This validates the SBs to be used as background estimators for the ROI. 

4. Results 

The single-energy spectra around the respective ROI are shown in Fig. 6 for all decay 
modes. All two-detector coincident events with decay mode optimized individual 
detector threshold and sum energy limit are shown in light gray; no other cuts are 
applied. The corresponding MC background model is shown in black for illustration. 
Events passing all cuts have exactly one entry in one of the 5 intervals (but two entries 
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Table 3. Summary of SBs after the complete sequence of cuts for each decay mode. 
Listed are the relative position A E of the four SBs compared to the ROI, the number 
of events Nmc in the MC background model and the number of events Ndata in the data 
set. Additionally the background expectations in the ROI and the observed events are 
shown after application of all cuts. The uncertainty of the MC expectations denote 
statistical uncertainties only. 


Region 

A E 

Nmc 

Ndata 

Region 

A E 

Nmc 

Ndata 

0+ s . -2+ 559.1 keV 

0 g + s . - 0+ 559.1 & 563.2 keV 

SB 1 

-7.5 keV 

2.5 ±0.1 

2 

SB 1 

— 12 keV 

8.0 ±0.1 

7 

SB 2 

+7.5 keV 

2.4 ±0.1 

0 

SB 2 

±12 keV 

7.5 ±0.1 

11 

SB 3 

— 15keV 

2.6 ±0.1 

3 

SB 3 

—24keV 

8.3 ±0.1 

7 

SB 4 

+15 keV 

2.4 ±0.1 

5 

SB 4 

±35 keV 

6.8 ±0.1 

9 

average SB 


2.5 ±0.1 

2.5 

average SB 


7.7 ±0.1 

8.5 

ROI 


2.5 ±0.1 

2 

ROI 


7.9 ±0.1 

5 

0+s. - 2+ branch 1: 559.1 & 657.0 keV 

0+ s . — 2\ branch 2: 1216.1 keV 

SB 1 

—8 keV 

8.5 ±0.1 

6 

SB 1 

— 19 keV 

0.52 ±0.02 

1 

SB 2 

+18 keV 

8.1 ±0.1 

5 

SB 2 

±10 keV 

0.39 ±0.02 

0 

SB 3 

— 16 keV 

8.7 ±0.1 

6 

SB 3 

—27keV 

0.58 ±0.02 

0 

SB 4 

+35 keV 

7.9 ±0.1 

12 

SB 4 

±47 keV 

0.32 ±0.02 

1 

average SB 


8.3 ±0.1 

7.25 

average SB 


0.45 ±0.01 

0.5 

ROI 


8.3 ±0.1 

6 

ROI 


0.40 ±0.02 

0 


in the single-energy spectrum). The ones in the ROI are shown in red and the ones of 
the SBs in blue. 

Frequentist 90 % confidence level and Bayesian 90 % credibility lower values were 
calculated for Ti/ 2 . The Frequentist values were based on a bounded profile likelihood 
test statistic [41] increased by 2.7 compared to the minimum. It was verified that this 
method has always sufficiency coverage. For the Bayesian credibility limit, a flat prior 
in TP g was assumed as well as a flat prior in the background level. 

In both approaches the same likelihood was used which is constructed for the inverse 
half-life Tyl. In case of two decay branches as for 2 u/3/3 0+ s — 2f, the likelihood is treated 
as two individual data sets with a common T~f. The expectation for the signal counts 
for a given decay branch (. k ) is: 

s k = In 2 • rj k ■ £ ■ Ty\, (3) 

with the decay branch specific efficiency ?%. The efficiency is the live-time weighted 
averaged detection efficiency of a 2 u/3/3 decay excited state event over all runs in the 
two array configurations of Gerda Phase I. The exposure £ is defined as the combined 
total isotopic exposure for 76 Ge in the data set for all detectors. The expected number 
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Figure 6 . Single-energy spectra around the respective ROI for all decay modes. 
Shown are all two-detector events for the optimized individual detector threshold and 
sum energy limit (light gray) and the corresponding background curves (black). The 
optimized cuts result in different two-detector spectra for each decay mode. Also shown 
are the ROI (shaded red) and SB region (shaded blue). Highlighted are events that 
are tagged as ROI (red) and SB (blue) after all cuts and that are used for the limit 
setting. Note that the histograms contain two entries per event and that one entry 
may lie outside the tagging region. 


of events in the ROI is: 

bk , 

t^k ^ Sfc 


(4) 


using the total background from all 4 SBs b^. 

The full expression of the likelihood is constructed with three terms: (1) a Poisson 
term describing the probability of signal plus background in the ROI, (2) a Poisson 
penalty term accounting for the uncertainty of the background level in the ROI and (3) 
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a Gaussian penalty term accounting for all systematic uncertainties condensed into the 
efficiency. The likelihood depends on T~^, the number of counts in the ROl (n k ), the 
number of counts in all SBs (m*,) and the efficiency expectation (e*,): 


C 


i ^ki Vk 


n 


r(w) <n,) --«,i 


r(« < ”** ) „-j 


1 1 ( tk-vk V 

1 . r A I 

L K)! J 


L (rn k )\ ' J 


eh fe v27r 


(5) 


For the profile likelihood we calculate —21og£. In the extraction of the posterior 
probability [42], flat priors for all fit parameters T~^, bk and 7jk were used since the prior 
information on the background and efficiency is included in the likelihood with penalty 
terms. 

Systematic uncertainties on the signal efficiency are estimated with MG simulations 
and combined into a single value cr £fe assuming no correlations. The sources of these 
uncertainties are the active volume and dead layer thicknesses of the detectors (5 %), the 
energy resolution and energy shift after cross talk correction (3%), the MC simulations 
(4%) and the uncertainty on the isotopic abundance (2.5%), where the numbers in 
parentheses give the resulting uncertainty on the signal efficiency. The combined relative 
systematic uncertainty on the efficiency is 7.5 % for all decay modes. The systematic 
uncertainty was investigated for the 0 +g — 0 ^ transition and assumed to be similar for 
all decay modes. 

A sensitivity study was performed using a toy MC under the assumption of no 
signal. The inverse half-life limit was calculated 10 4 times with randomly changing input 
parameters; ?% and and ??% and were each randomized according to a Poisson distribution. 
The expectation values of this distribution were taken from the background model 
prediction Nmc and 4Nmc, respectively (Table 3). e^, and is randomized with a Gaussian 
distribution with a mean of e^ and a width of a ek . The median of the 90% quantile 
inverse half-life limit distribution is taken as the sensitivity. 

Table 4 shows the input parameter for the statistical analysis (n k , irik and e k ) in 
the first 3 columns. The next two columns show the Frequentist 90% C.L. lower value 
for T 1/2 and the expected sensitivity T \/ 2 from the toy MC study. The last two columns 
show the Bayesian 90% credibility lower bound on Tu 2 and the respective sensitivity 

Td\ 1‘1- 


decay mode 2 u/3/3 0+ s — 2f: The transition to the excited 2f has one 559.1 keV de¬ 
excitation 7 ray and hence a small coincidence efficiency compared to the other excited 
state decay modes. After all cuts 2 events are observed in the ROI and 2.5 events are 
expected from the side bands. No signal is found and the Frequentist analysis yields a 
90% lower value on the half-life: T'A 2 > 1.6T0 23 yr. The sensitivity as defined above is 
1.3-10 23 yr. The Bayesian analysis yields lower credibility bounds on the half-life of TfY 2 
> 1.3-10 23 yr (90% C.l.) with a sensitivity of 1.2-10 23 yr. 
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Table 4. Summary of results for all decay modes. Shown from left to right are 
the input parameters for the likelihood: n - number of events within the ROI, m - 
number of events with in SB (1-4), e - detection efficiency. The following columns show 
the Frequentist lower half-life limit and the sensitivity which contains the true value 
in 90 % of the cases. The last two columns show the Bayesian lower limit for half-life 
values with more than 90 % probability and the respective sensitivity. 






Frequentist 90 % C.L. 

Bayesian 90% C.I. 

Decay mode 

n k 

rrik 


Ti /2 

T-L/2 

Ti/2 

Ti/2 




[%] 

o 

to 

CO 

. 1-3. 

o 

to 

CO 

o 

to 

CO 

o 

to 

CO 

. i-i . 

- 2 f 

2 

10 

0.389 

> 1.6 

> 1.3 

> 1.3 

> 1.2 

o+. - of 

5 

34 

0.919 

> 3.7 

> 1.9 

> 2.7 

> 1.8 

0 +s. — branch 1 

6 

29 

0.594 

> 1.7 

> 1.2 

> 1.4 

> 1.1 

0 +s. — 2 ^ branch 2 

0 

2 

0.092 

> 0.74 

> 0.64 

> 0.49 

> 0.46 

0 + s . — 2f combined 

- 

- 

- 

> 2.3 

> 1.4 

> 1.8 

> 1.3 


This lower half-life limit for the 2 u/3/3 0^ s — 2f transition is two orders of magnitude 
better than previous best limits. However, the current theoretical prediction are beyond 
the experimental reach. The most recent calculation [28] predicts a half-life longer by 
4 orders of magnitude. The lowest half-life prediction by the Hartree-Fock-Bogoliubov 
approach is still a factor of 6 above the current sensitivity. The largest half-life prediction 
by shell model calculations is even 7 orders of magnitude above the current limit. 

decay mode 2 u/3/3 0^ s — 0^: The transition to the excited 0^ state has two de¬ 

excitation 7 rays of 559.1 keV and 563.2 keV in the final state and the largest efficiency 
of the considered decay modes. It also has the lowest theoretically predicted half-life 
and is thus of special interest. After all cuts, 8.5 events are expected in the ROI and 5 
events are observed. No signal is found; a lower half-life value of 3.7-10 23 yr (90% C.L.) 
is set for the Frequentist analysis and 2.7-10 23 yr (90% C.I.) for the Bayesian analysis. 
The sensitivities are 1.9-10 23 yr and 1.8-10 23 yr, respectively. 

This lower half-life limit is 2.5 orders of magnitude better than previous results for 
the 2v(3f3 0+ s — 0^ transition. The new limit is well within the region of theoretical 
predictions. Bayes factors are calculated for testing the hypothesis of each NME model 
in Table 1 by taking the ratio B = p(Hi)/p(H 0 ) in which Hi is the NME model 
hypothesis with T^Jf del and H$ the hypothesis of only background. The models in 
Refs. [22, 23, 24, 25] have B <10 ~ 6 and are ruled out. The QRPA model [26] has 
B = 0.001 — 0.19 for T^ = (1.0 — 3.1)-10 23 yr, respectively. Recent calculations with 
RQRPA [11], IBM-2 [14] and SM [16] predict significantly longer half-lives. For RQRPA 
a range can be constrained: B = 0.005 for c/a = 1-00 (T 1 , 2 = 1.2-10 23 yr) compared to 
B = 0.45 for (]a = 1-26 (T ^ 2 = 5.8-10 23 yr). The IBM-2 and SM prediction are still 
above the current experimental reach. 
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decay mode 2 u/3/3 0+ s — 2f: The transition into the excited state has two 

de-excitation branches: Branch 1 (64% probability) with two 7 ray emissions of 

559.1 keV and 657.0 keV and branch 2 (36 % probability) with a single 7 ray emission of 

1216.1 keV. Branch 1 shows a significantly larger efficiency due to the higher branching 
ratio; on the other hand, the ROI for branch 2 is at higher energy resulting in lower 
background level. 7.25 events are expected and 6 events are observed for branch 1 
compared to 0.5 expected and no observed event in branch 2. No signal is found in 
either branch and a combined limit is calculated according to Eq. 5 with k = 1,2: 
The Frequentist 90% C.L. lower value is 2.3-10 23 yr. The Bayesian analysis yields a 
90% credibility lower bound of 1.8-10 23 yr. The corresponding half-life sensitivities are 
1.4-10 23 yr and 1.3-10 23 yr, respectively. 

For the 2 u/3/3 2^ s — 2f transition the lower half-life limit was improvement by 2 
orders of magnitude compared to previous limits. The theoretical predictions for this 
transition are > 10 28 yr and cannot be tested with the current sensitivity. 

5. Conclusions 

An analysis for 2 10/3/3 excited state transitions in 76 Ge with the Gerda Phase 1 data 
set has been performed for the three decay modes 0 + s — 2f , 0 + s — 0 ^~ and 0 + s — 2%. 
The analysis is performed without blinding, however the automated choice of cuts is 
expected to have reduced bias: All cut parameters are chosen such that the sensitivity 
calculated from MC simulations is maximized. No signal has been found and new half- 
life lower limits are set for all decay modes which are at least two orders of magnitude 
larger than those reported previously. Bayes factors are calculated for the predictions 
of the 0+ s — 0^ half-life with various nuclear models. Many old NME calculations could 
be ruled out. 

The analysis is based on the assumption that only one decay mode is realized at 
a time. This is valid in the present case for the non-observation of a signal. Hence, 
the analysis is performed on each decay mode completely independently. However, it 
should be noted that the results of the different decay modes are not decoupled since 
they proceed through the same levels. Particularly the 559.1 keV 7 line of the 0+ s — 2\ 
transition is part of all decay modes. The Frequentist lower half-life limits are larger 
than the sensitivity in all cases. A statistical background downward fluctuation in the 

559.1 keV region, as observed, has a similar influences on all limits. For additional 
information refer to a more detailed description of this analysis in Ref. [40]. 

In the future Phase II of the Gerda experiment it will be possible to increase the 
sensitivity further. The target mass of enr Ge detectors will be increased by a factor 
of two in form of relatively small BEGes detectors. The overall background level is 
expected to be reduced by an order of magnitude. 
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